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ABSTRACT 

Liquid  diffusion  coefficients  for  systems  approxi¬ 
mating  gas  chromatographic  sample  —  liquid  phase 
systems  were  determined  by  using  multiple  internal  re¬ 
flectance  infrared  spectroscopy.  An  expression  for  the 
diffusion  coefficient  as  a  function  of  the  change  of 
concentration  with  time  was  developed  from  FickTs 
second  lav;  and  JostTs  treatment  for  diffusion  into  a 
slab.  By  using  cell  covers  designed  for  this  project 
(for  FMIR  solid  sampling  cell)  and  by  turning  the 
spectrophotometer  (Infracord)  on  its  side,  it  was 
possible  to  measure  the  change  of  concentration  of 
solute  (fed  as  vapour  in  carrier  gas)  in  a  layer  of 
solvent  with  time  spectrophotometrically  and  to  thus 
determine  the  diffusion  coefficient  of  the  solute- 
solvent  system.  The  systems  used  in  the  project  in¬ 
cluded  those  of  different  solutes  in  fluorolube  and 
acetone  in  solvents  of  varying  viscosities. 

An  attempt  was  made  to  qualitatively  evaluate  the 
data  by  correlation  with  known  values  and  values  ex¬ 
pected  from  empirical  expressions.  The  results  indicate 
that  this  method  provides  qualitative  values  for 
diffusion  coefficients  with  main  errors  due  to  the 
preparation  of  the  solvent  layer. 
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1  INTRODUCTION 


The  research  described  in  this  paper  arose  from 
interest  in  two  areas  of  instrumental  analysis,  namely 
gas  chromatography  and  FMIR  (Frustrated  Multiple  Internal 
Reflectance)  infrared  spectroscopy.  Interest  in  the 
liquid  mass  transfer  term  in  gas  chromatography  and 
subsequent  reading  led  to  the  conclusion  that  there  is 
very  little  data  available  concerning  diffusion  co¬ 
efficients  in  chemical  systems  applicable  to  chrom¬ 
atography.  Research  has  been  done  in  this  field  but  due 
to  the  lack  of  numbers,  only  qualitative  results  have 
been  obtained.  It  seemed  possible  to  use  the  technique 
of  reflectance  spectroscopy  to  measure  diffusion  co¬ 
efficients  of  systems  similar  to  the  sample  —  liquid 
support  systems  of  gas  chromatography  and  perhaps  provide 
some  numbers  for  the  mass  transfer  term. 

1-1.  Liquid  Phase  Mass  Transfer  Term  in  Gas  Chromatography 

In  gas  chromatography  column  efficiency  is  most  often 
expressed  in  terms  of  HETP  (Height  Equivalent  to  a 
Theoretical  Plate) .  The  value  of  H  for  a  particular 
system  is  usually  considered  to  be  dependent  on  three 
factors:  diffusion  spreading  (longitudinal  spreading  of 
the  plug  of  solute  molecules  in  the  mobile  phase) ,  eddy 
diffusion  (caused  by  variations  in  flow  velocity  at 
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different  points  in  the  packing  caused  by  the  packing) 
and  the  resistance  to  mass  transfer  in  both  the  mobile  and 
stationary  phase.  Van  Deemter,  Zuiderweg,  and  Klinkenberg 
(22)  expressed  this  dependency  mathematically;  in  its 
simpler  form  this  relationship  is 

H  =  2  +  2  VD  ,  2k'  d  2  u 

P  g  f  o 

uQ  3 (  1+k) 2D1 


or  as  more  commonly  written  H 
where  - 


A  +  B/u  +•  C ,u  ,  0-5) 
o  1  o 


the  average  velocity  of  the  carrier  gas 

the  eddy  diffusion  term 
the  particle  diamenter 

a  factor  relating  to  the  irregularities 

in  column  packing 

gas  diffusion  term 

gas  diffusion  coefficient 

tortuosity  correction 

liquid  phase  mass  transfer  term 

liquid  diffusion  coefficient 

distribution  coefficient  function 
-  average  thickness  of  the  liquid  phase 
on  the  solid  support. 


The  value  of  jct  is  defined  by  the  distribution  coefficient 
K,  the  volume  fraction  occupied  by  the  gas  phase  F  ,  and 

—  o 

the  volume  fraction  occupied  by  liquid  phase  so  that  kT=r 
K  (Fx/Fg)  (20). 

The  mass  transfer  term  as  given  above  was  derived 
(21  )  without  any  special  regard  to  the  type  of  column  or 
application  of  solvent  in  the  column.  One  special  case, 
which  may  occur  in  packed  columns  and  which  is  generally 
assumed  in  capillary  columns  is  that  of  a  continuous  film. 

The  term  defined  by  Golay  for  columns  &3  )  is 

k^  u  r  ^ 

0  0 

6  (1+k)2  F2c2D1 

where  rQ  is  the  radius  of  the  capillary,  k  is  the  distri¬ 
bution  coefficient  terra,  F  is  a  constant  of  proportionality 
used  in  the  determination  of  effective  area,  and  c  is  the 
partition  coefficient. 

In  both  expressions  for  C1  it  is  seen  that  the  mass 
transfer  term  is  directly  proportional  to  the  square  of 
the  thickness  of  the  layer  of  liquid  phase  and  to  the 
velocity  of  the  carrier  gas;  and  it  is  inversely  pro¬ 
portional  to  the  diffusion  coefficient  in  the  liquid  phase. 

Work  has  been  done  in  separating  the  term  from  the 
general  rate  equation  and  finding  experimental  values  for  it. 


. 
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Barr  and  Sawyer  (1)  found  the  functional  dependence  of 

on  k  (or  kT ) ,  and  found  the  effect  of  different  carrier 

gases  on  the  value  of  the  k  terra.  Perrett  (19)  in  studying 

the  effect  of  column  repacking  on  the  values  of  C,  and  C 

1  g 

(where  is  the  term  corresponding  to  the  resistance  to  mass 
transfer  in  the  gas  phase  used  in  the  extended  Van  Deemter 
equation)  found  values  of  k  and  C-^  and  from  these  cal¬ 
culated  values  of  df2/D1  (Table  1).  To  calculate  the 
values  of  d^/D-^  he  assumed  the  Van  Deemter  term  depen¬ 
dence  of  C-^  on  k  as  f  (k)  =■  2k/3  ( 1  k)2.  He  states  that  this 
f(k)  seemed  to  give  reasonable  values  of  d  2/d  but  that 
the  lack  of  independent  diffusion  coefficient  data  provided 
only  qualitative  results.  Bohemen  and  Purnell  (3)  calcu¬ 
lated  values  of  C-^  using  various  particle  sizes  of  solid 
support  (Table  11)  and  finding  a  constant  value  of  d^2/D 
for  all  runs  concluded  that  for  fine  particles  d^  was  in¬ 
dependent  of  particle  size.  In  arriving  at  these  con¬ 
clusions  these  authors  assumed  the  value  of  the  diffusion 
coefficient  of  acetone  in  polyethylene  glycol  to  be 
10”?  cm^sec'^.  Evanoff  (7)  studied  the  effect  of  solvent 
and  temperature  on  diffusion  coefficients  of  some  organic 
nitro-compounds  using  polarography  and  viscosity  measurements. 

The  work  quoted  above  is  representative  of  what  has 
been  done  and  indicates  the  lack  of  independent  diffusion 
coefficient  data  for  these  systems.  For  this  reason  this 
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TABLE  1.  VALUES  OF  C1  AND  df2/D1  FOUND  BY  ELUTING  PROPANE, 
ISOBUTANE  AND  n-BUTANE  ON  A  COLUMN  OF  n-OCTADECANE  ON  SIL-O- 
CEL  C22  MESH  RANGE  100-120  BSS  USING  AND  N^  AS  CARRIER 

GASES  (19) 


Column 

Carrier 

gas 

X 

i — i 

O  c°\ 

O 

104, 

^4 

sec 

n-C 

4 

df2A 

c 

3 

>1  x  10' 
i-C 

4 

3 

,  sec 

n-C 

4 

1  A 

H 

2 

13*1 

l6tl 

1U1 

8.1*1 

13.5tl 

11.8+1 

N2 

16 

15 

12 

10 

12.6 

12.9 

I  B 

H2 

16 

14 

11 

10 

11.8 

11.8 

N2 

16 

15.5 

11 

10 

13.0 

11.8 

I  C 

H 

2 

15 

16 

11 

9.4 

13.4 

11.8 

13 

16 

12 

8.1 

13.4 

12.9 
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TABLE  11.  VALUES  OF  Q  FOUND  BY  ELUTING  ACETONE  AND  BENZENE 
ON  COLUMNS  WITH  VARYING  SUPPORT  SIZE  OF  POLYETHYLENE  GLYCOL3, 
ON  SIL-O-CEL  USING  H2  AND  N2  AS  CARRIER  GASES  (3) 


Sil-O-Cel  Support 

Cp  x 

10^,  sec 

Size,  dp,  cm 

Carrier  gas 

Acetone 

Benzene 

O.Q32 

H2 

2./, 

2.9 

N2 

2.4 

2.9 

0.016 

H2 

1.1 

1 . 6 

N2 

1.1 

1.6 

0.010 

H2 

1.2 

1.5 

N2 

1.2 

1.5 

do  percent  W/W 


t 
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work  was  begun  to  try  to  find  a  way  of  providing  numbers 
which  could  give  more  information  about  phase  thicknesses, 
about  C-j^  itself  and  its  functional  dependence  on  k. 

1-2 .  Diffusion 

Diffusion  is  a  process  which  leads  to  the  equali¬ 
zation  of  concentrations  within  a  single  phase.  The 
lav/s  of  diffusion  connect  the  rate  of  flow  of  the  dif¬ 
fusing  substances  with  the  concentration  gradient  re¬ 
sponsible  for  this  flow.  The  principle  rules  which 
govern  diffusion  are  FickTs  lav/s  (6).  FickTs  first  law 
states  that  the  flux  J  of  a  given  component  or  solute 
across  a  unit  area  fixed  in  the  medium  is: 

-D  i c/  c*x 

2  i 

where  D  is  the  diffusion  coefficient  in  units  cm  sec 
and  c/ ^  x  is  the  concentration  gradient.  Often  it 
is  not  possible  to  investigate  diffusion  coefficients 
under  conditions  of  a  constant  concentration  gradient; 
instead,  it  is  preferable  to  study  the  change  of  concen¬ 
tration  with  time  within  the  medium.  When  there  is 
diffusion  in  the  x-direction  only,  and  if  the  rate  of 
increase  in  the  amount  of  substance  within  a  volume 
element  bounded  by  two  parallel  planes  of  unit  area  situ- 
atex  at  x  and  x-f-dx  is  observed,  this  increase  is  equal 


* 


w  ■  ' 
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to  the  amount  diffusing  into  the  volume  element  minus  the 
amount  diffusing  out  of  it  per  unit  time,  that  is 

Jx  “  ^x-k}x'  Thus  the  rate  of  accumulation  in  the  volume 
element  is: 

JX  -  Jx+dx=  "D  (  ic/  ^x)x  ■  (_D  (  *°/  ^x)x+dx] 

—  -  (  ^  J/  ^>x)  dx . 

The  rate  of  concentration  increase  ^•c/'kt  is  the  accumu¬ 
lation  rate  per  unit  volume  ( J  —  J  )/dx,  i.e., 

X  X>pdX 

^c/  *t  f=  -  3j/  <*x. 

Substituting  for  J  from  FickTs  first  law  the  expression 

^  c/  <^t  =  D  V^c/  ^x^ 

is  obtained,  which  is  FickTs  second  law  of  diffusion  (  9) 
In  Fick’s  second  law  the  diffusion  coefficient  is  assumed 
constant . 

The  experimental  system  in  this  work  is  approximated  by  a 
slab  bounded  by  two  planes  perpendicular  to  the  x-axis 
(either  of  infinite  dimensions)  or  bounded  by  an  imper¬ 
meable  surface  of  arbitrary  shape  but  parallel  to  the 
x-axis.  In  this  bounded  system  the  initial  and  boundary 
conditions  are  only  dependent  on  x. 

For  diffusion  into  a  slab  Jost  (16  )  begins  the 
solution  of  FickT s  second  law  by  putting  concentration  as 
a  function  of  time  and  distance,  i.e., 
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c  -X(x)  T(t)  .  (i) 

This  makes  the  second  law 

X  T’  -  D  T  X”  (ii) 

where  the  primes  refer  to  the  differentiation  with  re¬ 
spect  to  the  corresponding  variables.  Thus  it  can  be 
written  in  the  form: 

T'/T  =  QX"/X.  (iii) 

These  can  only  be  equal  if  they  are  equal  to  the  same 

2  o 

constant  which  is  chosen  to  be  -  ^  D.  Thus,  Tt/T~  -  X  D 

and  Xn/X  ~  -A2  with  solutions  to  the  differential 
equations  being 

T  ~  exp  (  -  X2Dt)  (iv) 

.X  —  A  cosXx  -h  B  sinXx.  (v) 

If  the  solution  is  to  be  finite  for  all  values  of  t, 

2 

then  X  must  be  greater  than  zero. 

For  the  unbounded  (infinite)  linear  system,  A  and 
B  are  considered  to  be  functions  of  A  ,  and  expression 

(i)  is  integrated 

<*> 

c  \  A(  A)  cosA*  +  B(  A  )  sin  Ax  exp  (-A2Dt)  d'X  . 

J  (vi) 

© 

For  the  bounded  system  the  concentration  is  again 
considered  to  be  a  function  of  time  and  distance  and 
solutions  (iv)  and  (v)  are  again  obtained.  Where  in  the 
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infinite  case  A  was  restricted  only  by  the  condition 
2 

v>  0,  in  the  bounded  case  it  is  restricted  further  and 


only  allowed  positive  integral  values.  The  integral  for 
c  is  now  represented  by  an  infinite  sum  over  all  discrete 
values  of  X  instead  of  as  before,  by  an  integral  over 


the  continuously  varying  values  of  X  • 


Specific  examples  for  diffusion  out  of  and  into  a 
slab  are  treated  here  because  of  their  special  applica¬ 
bility  to  this  work. 

For  the  diffusion  out  of  a  slab  it  is  necessary  to 
set  up  the  proper  boundary  conditions  which  are: 
c  =  cQ  for  0<  x<"h  at  t  -  0, 
c  =•  0  for  x  ="  0  and  xrh  at  t^O 
(i.e.,  diffusion  from  both  sides  of  the  slab). 

2 

From  expression  c  ^  (A  sin  Ax  4-  B  cosXx)  exp  (-A  Dt) 
it  is  seen  that  the  boundary  conditions  are  fulfilled 
for  B=  0  and  7\  =  nTr/h. 

The  general  solution  is  now 


sin  nTT  x/h .  ( vii) 


The  arbitrary  constants  An  must  be  determined  in  such 
a  wav  as  to  conform  with  the  initial  conditions,  i.e., 


for  0<^x<^h  at  t=0 .  (viii) 


' 
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Jost  uses  Fourier's  theorem  and  arrives  at  the  solution 

to 

c  ~  /fCo/lT  L.  (l/2v+l)  sin  (2y+l)  It x  exp  -  f(2v+l^ 

\j-o  h 


v  -  0,  1,  2 ,  3  . .  . 


Dt 


(ix) 


.  When  the  solution  is  plotted  as  c  versus  x  for 
different  times,  it  is  seen  (17)  that  the  curves  are 
symmetrical  around  h/2,  which  is  to  say  that  there  is  no 
diffusion  across  the  midpoint  of  the  slab. 


For  the  problem  of  diffusion  into  a  slab  the  solu¬ 
tions  are  easily  derived  from  the  above  by  consideration 
of  the  new  boundary  conditions  i.e.,  c  —  cQ  at  x  -  0 
and  at  xrh,  for  t  greater  than  zero.  If  (ix)  is  a  solu¬ 
tion  of  FickTs  second  law,  the  negative  of  it  is  also  a 

solution.  If  c  is  a  solution  as  in  ( ix) ,  then  c  -  c  is 

o 

the  solution  of  the  negative  problem.  Thus  the  solution 
for  diffusion  into  a  slab  is  given  as 


c  —  c 


C  00 

1  -  44  2  1 

V  Vs£>  2vv 


1  sin  (2v41)TTx  exp 
1  h 


( 2v+l)ri)  Dt 


3 


If  c  were  to  be  plotted  versus  x  at  various  times, 
the  plot  would  resemble  that  shown  in  fig.  1,  and  be 
symmetrical  around  the  mid-point  of  the  slab  (i.e.,  at  h/2) 


J 


The  same  effect  could  also  be  caused  by  an  imper¬ 
meable  membrane  at  h/2.  Thus,  if  diffusion  through  only 


. 
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Figure  1.  Plot  of  integral  equation  (x)  for  diffusion 
into  a  plate,  the  concentration  at  the  boundaries  being 
kept  at  c 

o 
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one  face  of  the  slab  is  to  be  considered,  the  function  need 
only  be  taken  for  0<^x^h/2.  If  expression  (x)  is  approxi¬ 
mated  to  the  first  order  ( v  ^  0)  the  resulting  expression  is: 


c  ~  c 


,  0  x  (  h/2 .  (  xi) 


1  -  4/jj  sinTTx/h  exp  (7f/h)2  Dt) 

Since  in  this  v/ork  the  concentration  will  be  measured  at 
h/2,  the  expression  reduces  further  to 


c  ~  c 


1  -  4/fl-  sin  IT/ 2  exp  -^(  ir/h)  Dt) 


(xii) 


which  should  give  the  value  of  the  concentration  at  h/2 
at  any  time  t. 


To  make  the  expression  more  tractable,  it  is  differen¬ 
tiated  with  respect  to  t  and  the  following  expression  is 
obtained : 

dc/ dt  w  c  4TT D  e“  (ir/h)2Dt.  (xiii) 

2 

h 

To  remove  the  exponential  logarithms  are  taken,  obtaining 

2  2 

In  (dc/dt)  =  In  (K)  +  (-TT  /h  )  Dt  (xiv) 

which  reduces  to 

log  (dc/dt)  =  (log  K)  -  ^  2Pt  (xv) 

2.3  h2 

in  which  log  (dc/dt)  is  a  linear  function  of  t. 

Thus  if  plots  of  concentration  versus  time  are 
taken,  their  slopes  at  various  values  of  t  can  be 
measured;  and  if  the  base  ten  logarithms  of  these  slopes 


' 


Y~ 
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are  plotted  against  time,  the  result  should  be  a  straight 
line  from  whose  slope  the  diffusion  coefficient  can  be 
calculated. 

Diffusion  coefficients  have  been  calculated  from  ex¬ 
pressions  relating  diffusion  to  more  measureable  physical 
properties,  especially  viscosity.  The  first  equation 
which  had  any  success  was  the  Stokes-Einstein  equation  (10): 

D1^RT/6'TT  r±  ^  N 

where  r^  is  the  "radius"  of  the  solute  molecule,  N 
is  AvogadroTs  number,  ^  is  the  viscosity  in  centipoise, 

R  the  gas  constant,  and  T  the  absolute  temperature.  This 
equation,  however,  was  derived  for  particles  of  size 
greater  than  the  solvent  molecules,  i.e.,  Brownian  (colloid) 
particles,  fohen  the  diffusing  particles  are  almost  the 
same  size  as  the  solvent  molecules,  a  slippage  factor  is 
introduced  and  the  6  replaced  by  4.  Many  other  empirical 
expressions  have  been  derived  such  as  the  Wilke-Chang 
equation  (ll)  (23) 

D_  —  7.4  x  lO-8  _  (xM) *  T 

1  v0*6 

where  x  is  an  association  parameter  which  is  introduced 
to  define  the  effective  molecular  weight  of  the  solvent 
with  respect  to  the  diffusion  process  (see  Table  111), 

M  is  the  molecular  weight  of  the  solvent,  T  is  the  absolute 
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TABLE  111.  WILKE-CHANG  ASSOCIATION 
P ARAI-ISTEHS  FOR  SELECTED  SOLVENTS  (23) 


Solvent 

Association  Parameter 

Water 

2.6 

Methyl  Alcohol 

1.9 

Ethyl  Alcohol 

1.5 

Benzene 

1.0 

Ether 

1.0 

Heptane 

1.0 
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TABLE  IV.  ATOMIC  VOLUKESaFOR  COMPLEX  MOLECULES  (23) 


Bromine  27.0 

Carbon  14.8 

Chlorine  24.6 

Hydrogen  3 ^ 

Iodine  37.0 

Nitrogen,  double  bonded  15.6 

Nitrogen,  in  primary  amines  10. 5 

Nitrogen,  in  secondary  amines  12.0 

Oxygen  (except  as  noted  below)  7.4 

Oxygen,  in  methyl  esters  9.1 

Oxygen,  in  methyl  ethers  9.9 

Oxygen,  in  higher  esters  and  ethers  11.0 

Oxygen,  in  acids  12.0 

Sulphur  25.6 

For  three-membered  ring,  as  in  ethylene  oxide,  deduct  0.6 
For  four-membered  ring,  as  in  cyclobutane,  deduct  $.5 

For  f ive-membered  ring,  as  in  furan,  thiophene,  deduct  11.5 
For  pyridene,  deduct  13 

For  benzene  ring, deduct  15 

For  naphthalene  ring,  deduct  30 

For  anthracene  ring,  deduct  47.5 


a.  atomic  contributions  to  the  molal  V3lume  (at  normal 
boiling  point)  for  complex  molecules. 
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temperature  and  V  is  the  molal  volume  of  the  solute  at 
normal  boiling  point  measured  in  cc/gm  mole.  The  value 
of  V  in  the  equation  is  estimated  from  Table  IV. 

In  this  work  the  calculation  of  diffusion  coefficients 
relies  on  one  basic  assumption,  that  used  by  Golay  (13  )  : 
the  rate  of  transport  of  solute  molecules  across  the 
solvent  surface  is  considered  kinetically  faster  than  the 
diffusion  process  and  thus  has  no  effect  on  the  apparent 
value  of  the  diffusion  coefficient. 

1-3  Reflectance  Spectroscopy 

As  FickTs  second  law  expresses  the  diffusion  as  a 
function  of  dc/dt ,  it  remains  to  measure  this  quantity. 
Spectroscopy  is  useful  in  determining  concentrations 
through  the  linearity  of  the  Beer-Lambert  law,  which  re¬ 
lates  the  absorbance  of  a  species  to  its  concentration  by 
the  relationship  A  -  abc  ( A-absorbance ,  a-absorpt ivity , 
b-path  length,  c-concentration) .  Normal  transmission 
spectroscopy,  by  the  nature  of  its  sampling  methods,  would 
require  that  samples  in  a  diffusion  experiment  be  selected 
at  known  times  and  analyzed  individually  to  give  a  point-by¬ 
point  plot  of  concentration  versus  time.  Problems  such  as 
the  actual  selection  of  samples,  of  transfer,  and  of 
spectroscopic  sampling  would  seriously  affect  the  credibility 
of  the  results.  Ideally  the  diffusion  experiment  should 


' 
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take  place  in  the  sampling  area  of  the  spectrophotometer 
to  provide  both  immediate  and  continuous  analysis.  On  the 
basis  of  the  desired  conditions  for  the  experiment  it  was 
decided  that  reflectance  spectroscopy,  particularly 
Multiple  Internal  Reflectance,  would  be  the  most  useful 
technique  for  analysis.  The  reasons  for  its  usefulness  will  be 
discussed  after  a  note  on  the  theory  of  reflectance  spect¬ 
roscopy. 

Internal  reflection  may  be  compared  to  external  re¬ 
flection  as  follows:  in  external  reflection  the  incident 
beam  travels  through  a  medium  of  lower  refractive  index 
and  is  reflected  from  the  surface  of  a  medium  of  higher 
refractive  index.  In  internal  reflection  the  beam 
travelling  through  a  medium  of  higher  refractive  index 
is  reflected  from  the  surface  of  a  medium  of  lower  re¬ 
fractive  index.  Furthermore,  a  critical  angle  is  defined 
for  any  pair  of  media  as  the  angle  between  the  incident 
beam  and  the  normal  to  the  surface  of  the  reflecting 
medium  such  that  sine  ©c  -  nr/nd  (Qc  -  critical  angle, 

nr  -  refractive  index  of  the  rarer  medium,  n^  -  refractive 
index  of  the  denser  medium).  At  angles  less  than  ©c,  total 
internal  reflection  ceases;  at  angles  ©c  to  grazing,  total 


internal  reflection  occurs. 


. 
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With  both  internal  and  external  reflection  reten¬ 
tion  of  most  of  the  energy  of  the  incident  beam  is 
expected  after  multiple  reflections.  However,  Harrick 
0-4  )  found  for  external  reflections  from  a  silvered  surface, 
that  after  ten  reflections  there  was  a  retention  of  only 
3/5  of  the  total  power;  and  that  after  the  same  number  of 
reflections  from  a  semi-conductor,  he  found  a  retention 
of  3*7  x  10  initial  power.  With  internal  reflections 
he  found  that  for  a  silicon  surface  there  was  a  retention 
of  50  per  cent  after  more  than  100  reflections. 

A  phenomenon  of  total  internal  reflection .is  the 
penetration  of  the  incident  beam  and  subsequently  re¬ 
flected  beams  (in  multiple  internal  reflectance)  into  the 
rarer  medium  (that  with  lower  refractive  index)  (14). 

Maxwell Ts  theory  predicts  a  penetration  of  the  totally 
reflected  power  although  there  is  no  net  flow  or  pro¬ 
pagation  into  the  rarer  medium.  The  depth  of  penetration 
is  directly  proportional  to  the  wavelength  cf  the  light 
in  the  denser  medium  and  is  inversely  proportional  to  the 
ratio  n^n^  (see  fig. 2)  .  Because  of  the  penetration  of 
the  totally  reflected  power  across  the  surface,  it  is  con¬ 
venient  to  represent  the  path  of  the  radiation  as  shown  in 
Fig.  3.  It  must  be  remembered  that  this  is  not  the  actual 
path  of  the  radiation  since  there  is  no  propagation  into 
the  rarer  medium. 


.. 

' 


si 


snetrat ion 
Depth 
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Figure  2.  Penetration  deoth  into  the  rarer  medium  versus 
the  angle  of  incidence  at  a  wavelength  of  3*4  microns  ( 14) 


denser  medium 

TTTT^rrnw'fr  medium 


Figure  3.  Apparent  path  of  radiation  in  total  internal 
reflectance 


-21- 


Another  phenomenon  which  is  very  important  in  re¬ 
flectance  spectroscopy  is  the  dependence  of  the  refractive 
index  of  the  rarer  medium  (or  any  material)  on  the  wave¬ 
length  of  the  light  -  especially  at  an  absorption  band  in 
the  medium  (see  fig. 4) •  Thus  in  a  dielectric  plate  of 
refractive  index  n^,  when  the  beam  is  totally  internally 
relected  at  an  angle  near  the  critical  angle,  there  will 
be  a  certain  amount  of  penetration  into  the  medium, 
depending  on  the  angle  of  incidence.  If  the  rarer  medium 
absorbs  at  the  wavelength  of  the  incident  beam,  there  will 
be  an  absorption  of  the  energy  that  has  penetrated.  Since 
the  refractive  index  of  the  rarer  medium  changes  at  the 
absorption  band  (anomalous  dispersion)  (2  ) ,  the  ratio 
n^/n^  also  changes  and  the  critical  angle  will  change.  If 
the  original  angle  of  incidence  is  near  or  at  the  critical 
angle,  the  new  critical  angle  will  be  less  than  the 
original,  the  angle  of  incidence  will  be  greater  than  the 
new  critical  angle  and  total  reflection  will  cease.  If  the 
wavelength  is  scanned  past  the  absorption  band,  the  re¬ 
fractive  index  returns  almost  to  its  original  value  and 

/ 

total  reflection  is  resumed.  If  the  intensity  of  reflected 
light  is  followed  spectrophotometrically ,  a  reflectance 
spectrum  of  the  rarer  medium  is  obtained.  (25) • 

There  are  two  basic  forms  of  reflectance  spectroscopy 
in  the  infrared  region,  namely  Attenuated  Total  Reflectance 


a  ,  *i  ’Id 
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Refractive 

Index 

or 

Absorbance 


nPi 


solid  line:  refrac¬ 
tive  index  of  sample 

dotted  line :  absorp¬ 
tion  band  of  sample 


dashed  lines:  refrac 
tive  indices  of 
reflector  plates 


*■>  A  Increasing  wavelength  ' 

Figure  4-  Anomolous  dispersion  -  behaviour  of  refractive 
index  of  sample  at  an  absorption  band  (25) 


Abso rbance 


Increasing  wavelength 

Figure  5*  Internal  reflection  absorption  bands  resulting 
from  the  level  (n  n  n  )  of  the  plate  refractive  index 

P.»  p2,  p3 

relative  to  the  refractive  index  of  the  sample  (25) 
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( ATH)  and  (Frustrated)  Multiple  Internal  Reflectance 
(  (F)  MIR) .  The  first  work  in  ATR  was  done  by  Fahrenfort 
(8  )  in  1961  and  has  since  expanded  until  today  ATR  and 
FMIR  attachments  are  commercially  available  to  fit  in  most 
infrared  spectrophotometers  (single  beam  and  double  beam), 
including  the  Perkin-Elmer  IR-137  (Infracord) .  ATR  utilizes 
a  single  internal  reflection  in  a  prism  against  whose 
surface  the  sample  is  firmly  pressed.  FMIR  uses  a  long 
rectangular  prism  plate  with  up  to  25  reflections  against 
which  sample  may  be  placed'-  on  either  or  both  sides. 

(fig.  6) 

The  choice  of  prism  materials  is  one  of  the  most 
important  considerations  in  both  ATR  and  FMIR,  the  prism 
material  requiring  high  refractive  index,  high  purity, 
high  polish,  durability,  inertness,  and  of  course  no 
absorption  bands  in  the  region  being  used.  A  list  of 
some  of  the  most  common  materials  and  their  properties  is 
given  in  Table  V.  The  most  useful  material  in  FMIR  is 
KRS-5  (42  mole  %  Thallium  bromide  and  mole%  Thallium 
Iodide)  which  has  a  high  refractive  index  and  high  polish, 
but  the  disadvantages  of  solubility  in  base  and  toxicity. 
Silver  Chloride,  although  commonly  used  in  ATR  in  prisms 
and  pressed  plates,  is  soft,  pliable  and  easily  deformed 
under  pressure,  and  discolours  in  visible  light.  Germanium 
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Figure  6.  Schematic  diagram  of  FMIR  optical  system. 

A,  entrance  slit  of  spectrophotometer;  B,  source;  C  FMIR 
prism  plate;  D,  E,  F,  G,  adjustable  mirrors. 
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TABLE  V. 

PROPERTIES 

OF  PRISM  MATERIALS 

FOR  REFLECTANCE 

- 

SPECTROSCOPY  (  5~a) 

Material 

Refractive 

Index 

Transmission 
Range,  microns 

Comments 

NaCl 

1.519 

1-16 

low  n ,  brittle , 
use  for  organic 
liquids 

CsBr 

1.667 

1-36 

long  range 
transmission 

AgCl 

1.997 

2-22 

soft,  discolours 
easily 

KRS-6 

2.192 

2-20 

intermediate  n, 
tough 

Irtran-2 

2.245 

2-13 

very  durable 

KRS-5 

2.380 

2-35 

most  useful, 
toxic,  soluble  in 
base  . 

Ge 

.4.017 

2.5.-  15 

highest  index 

. 


. 
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is  especially  good  for  samples  loaded  with  free  Carbon. 

KRS-6  (40  mole  %  Thallium  Bromide  and  60  mole  %  Thallium 
Chloride)  is  extremely  useful  for  a  long  wavelength  range. 
Irtran-2  has  good  physical  characteristics,  i.e.,  high 
refractive  index,  inertness,  opaque  up  to  125°C,  but 
gas  high  reflection  losses  and  is  only  useful  in  a  short 
range  -  up  to  eleven  microns. 

Since  ATR  deals  with  a  single  reflection,  large 
penetration  is  desired  and  the  prism  material  refractive 
index  should  be  fairly  low.  This  is  shown  as  n  in  fig.  4. 
Unfortunately,  use  of  materials  with  a  lower  refractive 
index  (such  as  AgCl)  results  in  some  distortion  of  the 
reflectance  band  at  higher  wavelengths,  (fig. 5)  This 
distortion  is  not  found  in  FMIR  because  of  the  multiple 
reflections  which  then  require  less  penetration  per  re¬ 
flection,  allowing  a  prism  material  of  higher  refractive 
index  to  be  used.  Normally  a  material  with  refractive 

index  of  level  n  in  fig.  4. is  used. 

P3 

Spectral  distortion  can  also  occur  when  the  angle  of 

incidence  on  the  reflecting  surface  approaches  or  exceeds 

the  critical  angle  whose  sine  is  expressed  as  n  /n,  (or 

r  d 

in  terms  of  prism  and  sample  n  /n  ) .  Thus  for  a  given 

r  s  p 

prism  material  and  fixed  angle  of  incidence  caution  must 
be  used  in  the  choice  of  sample  material  and  its  associated 


i 
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refractive  index.  At  a  fixed  angle  of  incidence  for  a 
given  crystal  there  is  a  critical  index  above  which  total 
reflection  ceases.  That  is  to  say,  the  critical  index 
occurs  when  the  angle  of  incidence  equals  the  critical 
angle  of  the  system.  To  reduce  spectral  distortion,  a 
safe  index  is  defined  as  being  the  index  of  refraction  0.2 
units  below  the  critical  index  (24).  For  example  with  KRS-5 
(refractive  index  2.4)  the  critical  index  for  the  rarer 
medium  at  angle  of  incidence  30°  equals  2.4  x  sin  30°^  1.2. 
The  safe  index  is  then  placed  at  1.0.  Values  for  critical 
and  safe  indices  are  given  for  KRS-5,  AgCl  and  Ge  for 
various  angles  of  incidence  in  Table  VI. 

Reflection  spectra  are  not  identical  to  transmission 
spectra  but  are  similar  and  provide  valuable  spectroscooic 
information  of  otherwise  intractable  samples  such  as  ad¬ 
hesives,  rubbers,  fabrics  and  crystals  (5-b,  26).  Sampling  is 
simple,  the  only  criterion  being  the  availability  of  a 
smooth  surface,  and  the  results  are  independent  of  sampling 
thickne  ss . 

Since  in  reflectance  spectroscopy,  particularly  in 
FMIR ,  only  one  surface  of  the  sample  need  be  in  contact 
with  the  prism-plate  to  provide  a  good  spectrum,  it  is 
possible  to  leave  the  other  surface  open  to  the  air  or  to 
any  cell  applicable  to  the  desired  experiment.  This 
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TABLE  VI. 


CRITICAL  AND  SAFE  INDICES  FOR  KR3-5,  AgCl, 


AND  GERMANIUM  (24) 


Material 

Angle  of 
Incidence 

Critical 

Index 

Safe 

Index 

KRS-5 

0 

0 

co 

1.25 

1.05 

(n  =  2  .4) 

-r- 

Vn 

O 

! 

1 

1 

1.7 

1.5 

60° 

2.0 

1.3 

AgCl 

O 

O 

CO 

1.0 

o.s 

(n  -  2  .0) 

45° 

1.4 

1.2 

60° 

1.7 

1.5 

Germanium 

O 

O 

CO 

2.0 

1.3 

(n  —  4) 

45° 

2.8 

2.6 

-  “60° 

3.5 

3.3 
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property  makes  it  practicable  to  the  diffusion  experiment. 

In  this  experiment  a  film  of  constant  thickness  of  "solvent” 
can  be  placed  on  the  surface  of  the  prism-plate  and  then 
covered  with  a  cell  system  which  will  provide  a  constant 
flow  of  "solute".  If  the  film  is  of  considerable  thick¬ 
ness  relative  to  the  penetration  depth  of  the  beam  power, 
then  as  diffusion  takes  place,  the  change  of  concentration 
at  the  plate  surface  (or  h/2  of  the  theoretical  slab)  can 
be  measured  continuously  by  measuring  the  change  in  ab¬ 
sorption  at  a  fixed  wavelength.  By  choosing  a  viscous 
solvent  such  as  found  in  liquid  phases  in  gas.  chromatography , 
and  by  feeding  the  solute  as  a  vapour  in  a  carrier  gas,  a 
gas  chromatographic  situation  is  approximated.  The  change  of 
concentration  with  time  is  measured  and  from  this  the 
diffusion  coefficient  can  be  calculated,  and  information 
allowing  elucidation  of  the  liquid  phase  mass  transfer  term 
can  be  obtained. 

For  maximum  workability  of  this  method,  it  remains  to 
be  seen  whether  absorbance  values  found  by  Fl'I R  need  to  be 
transformed  into  concentration  values,  or  if  the  absorbance 
plots  can  be  used  directly  in  evaluating  the  diffusion  co¬ 
efficient.  Since  the  spectrophotometer  used  in  this  work 
records  on . a  linear  percent  transmittance  scale,  it  is 
preferable  to  replot  the  original  spectra  on  a  linear  ab¬ 
sorbance  scale.  Assuming  BeerTs  law  (A  —  abc)  to  oe  valid, 


" 
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the  slope  of  the  plot  at  any  time  t  can  be  expressed  as 
A  A/^t  ~  abA<f/£t,  where,  for  the  system  ab  remains 
constant.  Taking  the  logarithms  of  these  various  slopes 
the  expressions  become  log  ab  •+  log  (  A  c/  .At)  .  If  these 
values  are  plotted  and  the  slope  of  the  resulting  straight 
line  taken,  the  expression  for  the  slope  will  be 
A  log  (>6c/At)/At,  which  in  the  limit  would  be  Alog(  dc/dt )  /At . 
This  in  turn  is  the  slope  of  the  straight  line  expected  from 
equation  (xv)  page  '3  .  Thus ,  no  transformation  of  values 
from  absorbance  to  concentrations  need  be  made  and  the 
diffusion  coefficients  can  be  calculated  almost  directly 
from  the  spectra. 


11  EXPERIMENTAL 


11-1  Reagents 

The  reagents  used  in  the  experimental  work  as  solvents 
(layer  materials)  were:  Fluorolube  S-30  (Hooker  Electro¬ 
chemical  Company) ,  Carbowax  600,  Carbowax  1500,  Dinonyl- 
phthalate  (DNP) ,  D-C  Silicone  grease,  and  Ucon  50-HB  2^0X 
(all  from  F  &  M  Scientific  Company).  The  samples  chosen 
were  1-pentene  (Fisher  technical  grade) ,  2-methyl  butane 
(isopentane  Phillips  66  commercial)  diethyl  ether  (Malinck- 
rodt  anhydrous)  and  wash  acetone.  The  BeerTs  lav.'  plots 
(figs.  13  &  14)  were  made  using  dodecane  (Mathescn,  Coleman 
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and  Bell)  . 

11-2  Design  of  Apparatus 

A  Perkin-Elmer  IR- 137-Infracord-  with  a  Wilks  Model 
12  MIR  attachment  and  KRS-5  plates  was  used.  The  Infra¬ 
cord  was  equipped  with  an  outlet  to  an  external  recorder 
(Sargent  Model  SR  variable  speed) .  By  use  of  a  10  milli- 
volt  plug  the  scale  was  expanded  so  that  0-100  per  cent 
transmittance  used  full  recorder  scale. 

In  order  to  use  thin  films  of  constant  thickness 
in  FMIR  it  was  necessary  to  have  the  plates  in  a 
horizontal  position.  Since  in  normal  use  the  plates  are 
placed  vertically,  it  was  decided  to  turn  the  Infracord 
on  its  side.  A  dexion  support  was  bolted  to  the 
Infracord  as  shown  in  fig.  7  and  the  spectrophotometer 
was  turned,  braced  with  wood  and  wired  to  the  wall. 

The  dexion  allowed  the  controls  to  be  manipulated.  When 
operated  in  this  position,  the  Infracord  did  not  show 
any  serious  effects  with  respect  to  scan,  gain,  wavelength 
shift,  but  did  show  some  heating  effects  due  to  the 
alignment  of  the  source,  and  some  misalignment  of  the 
mirrors  causing  the  beam  to  be  shifted  slightly.  The 
mirrors  were  adjusted  and  to  counteract  the  heating 
effects,  the  machine  was  turned  off  and  allowed  to  cool 
every  two  or  three  hours . 
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Figure  7-  Schematic  diagram  of  the  spectrophotometer 
as  supported  on  its  side. 

a.  Tabs  are  welded  to  the  dexion  at  intervals  determined  by 
the  bolts  in  the  ledge  on  the  spectrophotometer.  Holes  are 
bored  in  these  tabs  and  the  dexion  is  bolted  to  the  spectro¬ 
photometer  through  the  tabs. 


4-3  rH 


Two  cell  covers,  CC-A  and  CC-B  (see  figs.  3,  9,  10) 
were  designed  for  the  FMIR  solid  sampling  cell  to  provide 
for  the  analysis  of  a  layer  of  constant  thickness  and  for 
passage  of  gas  over  the  layer.  CC-A  consisted  of  two 
pieces  of  stainless  steel  which  supported  the  prism- 
plate  on  both  sides  but  left  about  1.4  cm.  of  open  surface 
between  them.  Enough  space  was  left  between  the  bottom  of 
the  cover  and  the  top  of  the  plate  to  allow  insertion  of 
spacers.  These  spacers  were  constructed  from  shim  materials 
of  stainless  steel  and  brass  of  thicknesses  0.005cm.  -  0.023  cm 
(127  microns  -  5$4  microns) .  When  the  cell  was  assembled, 
the  plate  was  placed  on  the  solid-cell  base,  the  spacer  was 
placed  on  the  plate  and  the  cell  covers  in  turn  were  placed 
on  the  spacer,  leaving  an  area  3«4  cm.  x  1.2  cm.  of  plate 
open  to  the  surface.  The  liquid  could  be  spread  on  the 
open  surface  and  with  a  teflon  spreader  its  thickness  regu¬ 
lated  to  that  of  the  spacer.  In  all  runs  a  layer  of  KEL-F 
#90  (halof luorocarbon)  grease  was  smeared  on  the  bottom  of 
the  spacer  to  give  a  better  seal  between  the  spacer  and  the 
plate  without  interfering  with  the  absorption  measurements. 

CC-B  was  constructed  from  a  single  piece  of  steel  4.5  cm. 
x  3.7  cm.  x  1  cm.  with  a  cavity  1.2  cm.  x  3.4  cm.  x  0.#  cm. 
designed  to  fit  exactly  over  the  part  of  the  plate  exposed  by 

Two  holes  were  drilled  from  the  top  of  the  cell 


the  spacer. 


C  -  spacer,  D  -  FMIR  prism  plate.  All  measurements  are 
given  in  cm. 

a.  Holes  at  four  corners  all  drilled  as  shov/n  in  diagram. 
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Figure  9.  Schematic  Diagram  cell  cover  CC-B  exterior 
and  bottom  view:  A  -  cavity  0.#  cm  deep.  All  dimensions 
in  cm.  The  holes  are  drilled  as  in  figure 


—  - > 


0 

0 

11  — — ^ 

J 

u 

crt\ 


-> 


\ 

f  T 

i-VtY, 

l  V 

l/'  0>7 

f 


C^\ 


< - 3  i  - - -> 

O-i 


spot  welded  to 
cavity  roof 


Figure  10.  Spacer,  top  view,  and  baffle  system 


for  CC-B  cavity. 


-37- 


into  the  cavity  and  these  were  fitted  with  0.2  cm.  stainless 
steel  tuoing  sections  which  extended  15  cm.  beyond  the  cell 
and  parallel  to  it.  A  small  piece  of  metal  was  spot  welded 
to  the  roof  of  the  cavity  and  this  was  berfc  so  that  it  ex¬ 
tended  across  the  cavity  about  0.3  cm.  from  its  roof.  This 
"baffle"  prevented  disturbance  of  the  liquid  layer  by  a 
direct  jet  of  the  incoming  gas  (see  fig.  10).  This  cell 
cover  was  placed  directly  on  top  of  the  spacer  and  a 
second  layer  of  KEL-F  was  applied  between  the  spacer  and 
the  cover  to  give  a  good  seal. 

A  system  utilizing  two  four-way  stopcocks  was  designed 
to  feed  sample  to  the  cell  (CC-B) ,  i.e.,  to  provide  flow  of 
pure  carrier  gas  and  carrier  gas  plus  sample  vapour  (  see 
fig.  11)  .  The  carrier  gas  was  bubbled  through  a  fritted 
glass  bubbler  and  the  sample  saturator  was  immersed  in  a 
constant  temperature  bath  to  provide  constant  sample  vapour 
pressure-  This  system  kept  a  continuous  flow  through  the 
sample  preventing  pressure  build-up  in  the  flask  and  re¬ 
sultant  "explosion"  of  sample  vapour  in  the  system  when  it 
was  changed  from  carrier  gas  to  sample  flow. 

The  system  from  the  bubbler  to  the  cell  was  made 
basically  of  0.4  cm.  glass  tubing  with  joints  of  tygon. 

The  steel  tubing  of  the  cell  was  inserted  into  the  glass 
tubing  and  the  joint  sealed  with  black  wax.  The  cell  was 
attached  to  and  dismantled  from  the  system  at  the  first 
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Figure  11.  Sampling  system.  A,  four-way  stopcocks;  B,  Helium 
inlet;  C,  black  wax  seal;  D,  cell  withCC-B;  E,  4mm  glass  tubing; 

*  x 

F,  trap;  G,  bubbling  jar  (saturator);  H,  sample;  1,  fritted 
glass  bubbler;  J,  constant  temperature  water  bath,  K,  outlet  for 
waste  gas. 

The  solid  arrows  correspond  to  pure  carrier  gas  flow  through  the 
system:  the  dashed  arrows  correspond  to  the  sample  flow  through 
the  system. 


. 
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tygon  connection  so  as  not  to  disturb  the  wax  seal.  The 
sampling  system  was  positioned  near  to  the  cell  to  minimize 
the  amount  of  tubing  required  to  reach  the  cell  and  thus 
decrease  the  dead  volume  of  the  cell  and  the  time  required 
to  reach  saturation  of  sample  vapour  above  the  layer. 

11-3  Verification  of  Beer’s  Law  for  FMIR 

Before  proceeding  with  the  actual  measurement  of 
diffusion  coefficients  it  seemed  necessary  to  verify 
BeerTs  lav/  for  this  particular  system;  studies  of  ab¬ 
sorption  at  various  concentrations  were  made  using  CC-A. 
Solutions  of  dodecane  (vapour  pressure  -  less  than  1  mm. 
at  room  temperature)  ( 4  )  in  fluorolube  at  known  concentra¬ 
tions  were  prepared.  The  cell  was  assembled  with  254  micron 
spacer  (KEL-F  #90  grease  seal)  and  a  smooth  layer  of  the 
solution  was  applied  using  the  teflon  spreader.  The  cell 
was  inserted  in  the  spectrophotometer  sampling  area  and 
the  mirrors  of  the  FMIR  attachment  were  adjusted  to  focus 
maximum  beam  energy  on  the  entrance  slit  of  the  Infracord  — 
it  was  found  necessary  to  do  this  frequently  as  the  slightest 
movement  of  the  mirrors  could  cause  a  serious  energy  loss. 

The  wavelength  was  scanned  from  3-0  to  4.0  microns  (see  fig. 
12).  Since  fluorolube  does  not  absorb  in  this  region,  it 
was  not  necessary  to  use  fluorolube  on  the  reference  plate. 
The  resulting  peak  at  about  3*4  microns  (alkane  C-H 
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Figure  12.  Spectrum  analysis  for  calculation  of  Beer’s  Lav/ 
plot.  Spectrum  of  0 . 3 5M  solution  of  dodecane  in  fluorolube. 
A,  spectrum  as  it  appeared;  B,  spectrum  as  analyzed  for  peak 

height,  A  A. 
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stretching)  was  analyzed  by  measuring  the  peak  height  at 
the  point  of  maximum  absorption.  The  transmittance  units 
from  the  spectrum  were  transformed  into  absorbance  units 
by  Beer’s  law  A  -  log  l/T  and  the  resultant  values  of  A  A 
were  plotted  against  concentration.  For  data  and  plot 
refer  to  Table  Vll  and  fig.  13. 

The  dependence  of  absorption  on  the  path  length  was 
investigated  by  using  a  solution  of  known  concentration 
(0.26M)  of  dodecane  in  fluorolube  and  by  varying  the  area 
of  the  plate  covered  by  this  solution.  The  cell  was 
assembled  with  CC-A  and  the  254  micron  spacer.  The  area 
to  be  covered  with  the  solution  was  separated  from  the 
rest  of  the  plate  with  a  thick  smeab  of  KEL-F  across  the 
plate.  A  layer  of  the  solution  was  placed  on  one  side  of 
the  grease  smear  and  a  layer  of  pure  fluorolube  was  placed 
on  the  other.  The  cell  was  immediately  placed  in  the 
sampling  area  and  the  wavelength  scanned  -  speed  was 
necessary  in  order  to  avoid  mixing  of  solutions  across  the 
grease  barrier.  This  procedure  was  repeated  for  several 
coverages,  i.e.,  7 mm,  10mm,  21mm,  25  mm,  and  34mm.  The 
spectra  were  analyzed  in  the  same  manner  as  before  and  the 
resultant  absorption  data  were  plotted. against  the  effective 
coverage.  The  results  are  shown  in  Table  Vlll  and  fig.  14  . 

o 

11-4  Measurement  of  Diffusion  Coefficients 


For  the  measurements  of  diffusion  coefficients  the 
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TABLE  Vll.  ABSORPTION  AT  3-4  MICRONS  OF  SOLUTIONS 
OF  DODECANE  IN  FLUORCLUBE  AT  VARIOUS  CONCENTRATIONS 


Solution 

Concentration 

Ak 

A 

0.09  M 

0.026 

B 

0.13  M 

0.035 

C 

0.15  M  ■ 

0.036 

D 

0.17  M 

0.046 

E 

0.20  M 

0.051 

F 

0.22  M 

0.049 

G 

0.24  M 

0.060 

H 

0.26  M 

0.059 

I 

0.31  M 

0.075 

J 

0.35  M 

0.084 

K 

0.40  M 

0.090 

L 

0.44  M 

0.103 

0.12 
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Concentration  of  dodecane  in  fluorolube,  Molar 

Figure  13 .  BeerTs  Law  plot,  Absorbance  versus  concentration, 
for  solutions  of  dodecane  in  fluorolube  (Data  from  Table  Vlll) 
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TABLE  Vlll.  ABSORPTION  AT  3.4  MICRONS  OF  0.26  M 
SOLUTION  OF  DODECANE  IN  FLUOROLUBE  AT  VARIOUS 
COVERAGES  OF  FMIR  PLATE 


Coverage ,  mm 


7 

10 

21 

25 

34  (full  plate) 


0.013 

0.020 

0.035 

0.04^ 

0.056 


Plate  coverage,  mm 


1 

. 
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solute  solvent  systems  were  chosen  on  the  basis  of*  favor¬ 
able  physical  and  spectroscopic  properties.  The  solute 
samples  were  chosen  on  the  basis  of  vapour  pressure  -  of 

the  order  of  400  -  600  mmH  .  The  samples  chosen  were 

§ 

1-pentene  (vapour  pressure:  100  mm  at  -23.0°C,  400mm  at 
,  12 . $  ,  and  760mrn  at  30.1°)  f  2-methylbutane  (v.p.  :  100mm 

at  -20.2°,  400mm  at  10.5°,  and  760mm  at  27. g°)  and 
diethyl  ether  (v.p..:  100  mm  at  -11.5°,  400mm  at  17.9°  and 
760mm  at  34*6°)  (  4  ) •  The  liquid  layer  (solvent)  was 
chosen  to  be  fluorolube  because  of  the  lack  of  absorption 
bands  in  the  region  3-7  microns. 

The  wavelength  of  maximum  absorption  for  alkane  C-H 
stretching  (about  3.4  microns)  was  determined  by  scanning 
the  spectrum  of  pure  solute  and  the  spectrophotometer 
wavelength  drum  was  locked  in  that  position.  The  sample 
was  placed  in  the  gas  bubbling  jar  (saturator)  at  a  level 
well  above  the  fritted  glass  bubbler.  The  bubbling  jar 
was  placed  in  the  constant  temperature  water  bath  and  the 
thermostat  on  the  bath  was  set  for  20°C .  The  FMIR  cell 
was  assembled  with  CC-A  and  using  spacer  of  desired  thick¬ 
ness  the  layer  of  fluorolube  was  prepared.  After  removing 
CC-A,  an  even  layer  of  KEL-F  was  spread  on  the  bottom  of 
CC-B  and  this  cell  cover  was  then  placed  on  the  spacer  and 
plate.  The  cell  was  inserted  in  the  spectrophotometer  and 


. 
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the  inlet  and  exit  tubes  were  connected  to  the  sampling 
system  at  the  tygon  joints.  The  carrier  gas,  Helium,  was 
allowed  to  flow  through  the  system  with  flow  rate  0.55  ml/ 
sec.  With  pure  carrier  gas  in  the  system  the  recorder  was 
turned  on  and  the  time  scan  begun  to  establish  a  base-line. 
After  a  short  time  both  stopcocks  were  turned  simultaneously 
to  start  sample  flow  and  this  point  was  recorded  on  the 
chart  paper  as  zero  time.  Recording  speed  for  most  runs 
was  set  at  TtfastTT  or  five  inches  per  minute.  When  equili¬ 
brium  concentration  was  reached  in  the  layer  (when  there 
was  no  further  change  in  absorbance)  the  stopcocks  were 
again  turned  and  the  system  cleared  of  the  solute.  The 
flow  through  the  cell  was  then  cut  off,  the  cell  dismantled, 
cleaned  with  chloroform,  and  the  run  repeated.  For  each 
of  the  three  samples  forty  runs  were  made,  ten  for  each 
spacer  (four  spacers  of  thicknesses  127yu,  254^*, 
and  5$Aya  were  used)  . 

To  examine  the  effect  of  solvent  (liquid  layer)  on 
the  diffusion  coefficients,  experiments  were  carried  out 
using  D-C  Silicone  grease  (F  &  M  Scientific  Co.),  Carbowax 
600,  Carbowax  1500,  Dinonylphthalate ,  and  Ucon  -50-HB  2^0x- 
all  F  &  M  Scientific  Co.  products).  Infrared  spectra  of 
these  liquids  were  taken  and  compared  to  find  a  common 
region  of  no  absorption.  With  the  exception  of  dinonylph¬ 
thalate  these  liquids  showed  no  absorption  in  the  region 
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4  -  6  microns  (2500cm'1  -  1650cm'1).  Dinonylphthalate  had 
an  absorption  band  (C^O  stretching)  at  about  1730cm"1.  It 
was  decided  to  use  acetone  as  the  sample  and  its  band  at 
about  1750cm  1  to  monitor  the  change  in  concentration.  In 
the  acetone -dinonylphthalate  system,  the  slight  difference 
in  position  of  peak  maxima  allowed  the  chanme  of  concen¬ 
tration  to  be  followed  although  the  energy  of  the  beam  was 
seriously  reduced  by  the  absorption  of  DNP.  Using  only  the 
127  micron  spacer  for  the  very  viscous  solvents  -  Carbowax 
1500,  D-C  silicone  grease  -  and  both  127  and  3^1  micron 
spacers  for  DNP,  Ucon  and  Carbowax  600  (because  of  their 
lower  viscosity)  the  procedure  mentioned  above  was  re¬ 
peated  with  the  wavelength  selector  locked  at  5.75  microns, 
and  the  resulting  spectra  were  analyzed. 

11-5  Comparison  of  the  Rate  of  Sample  Delivery  to  the 
Surface  with  Rate  of  Diffusion  into  the  Liquid 

The  rate  of  delivery  of  sample  to  the  surface  of  the 
layer  must  be  large  compared  to  the  rate  of  diffusion  into 
the  liquid  for  the  diffusion  data  obtained  to  be  valid, 
i.e.,  concentration  of  solute  above  the  plate  must  be 
constant.  This  was  tested  by  measuring  the  amount  of 
solute  present  in  the  film  at  equilibrium  concentration 
as  a  fraction  of  the  solute  that  passed  through  the  cell 
cavity  during  the  time  taKen  for  the  layer  to  achieve 


. 


> 
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equilibrium  concentration. 

for  the  samples  of  higher  vapour  pressure,  1-pentene, 
2-methylbutane ,  diethyl  ether,  1-pentene  was  used  in 
showing  that  the  rate  of  depletion  of  sample  vapour 
from  the  cavity  is  in  fact  small  enough  to  Justify  the 
assumption  that  the  amount  of  sample  above  the  laver  re¬ 
mained  constant. 

Values  of  vapour  pressure  of  1-pentene  were  plotted 
(fig.  15)  at  several  temperatures  and  the  value  of  vapour 
pressure  at  20  c  was  obtained  from  the  resulting  curve  - 
550  mm  Hg. 

The  cell  was  assembled  with  CC-A  and  the  254  micron 
spacer;  a  layer  of  fluorolube  was  applied  to  the  plate. 

'With  1-pentene  in  the  bubbling  jar,  the  water  bath  temper¬ 
ature  set  at  20°c ,  the  flow  rate  at  O.34  ml/sec,  and  the 
wavelength  selector  of  the  spectrophotometer  locked  at  3*4 
microns,  the  cell  cover  CC-B  was  placed  on  plate  and  the 
procedure  for  measuring,  diffusion  coefficients  was  followed. 
The  point  corresponding  to  the  time  at  which  equilibrium  con¬ 
centration  in  the  layer  was  reached  was  marked  on  the  re¬ 
sulting  spectrum,  i.e.,  53  sec.  The  A  A  at  this  point  was 
calculated  and  from  the  BeerTs  law  plot  for  dodecane  in 
fluorolube  the  equilibrium  concentration  was  estimated. 

The  value  of  A  A,  0.05  absorption  units,  corresponded  to 
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concentration  value  of  0.2M.  The  amount  of  1-petene  in 
the  layer  was  calculated  to  be  1.2  x  3.4  x  0.0254  x 
0.2/1000  ^  2  x  10  5  moles .  The  time  required  to  deliver 
this  amount  into  the  sample  cavity  was  calculated.  From 
the  ideal  gas  law  PV  -  nRT ,  where  P  =  O.73  atm,  n  -  2  x 
10 moles,  R  —  0.0£  1  atm/deg  mole,  T  -  293°A,  and  V- 
•  (0 .34/1000)  t ,  the  time  was  calculated  to  be  two  seconds. 

In  any  run  there  will  be  a  time  required  for  the 
system  to  reach  an  equilibrium  concentration  of  sample  and 
since  this  is  related  to  the  dead  volume  of  the  cell 
(volume  of  delivery  system  from  stopcock  to  cavity)  this 
time  can  be  called  dead-time.  The  value  of  this  time  can 
be  estimated  by  following  the  procedure  for  measuring 
diffusion  coefficients  using  a  very  thin  layer  of  solvent 
(5-10  microns).  The  time  from  zero  to  time  of  initial 
change  as  shown  on  the  spectrum  is  designated  to  be  dead¬ 
time.  Using  1-pentene  and  a  very  thin  layer  of  fluorolube 
the  dead-time  for  the  system  is  seen  to  be  about  15  seconds 
(with  flow  rate  0.34  ml/sec) .  Thus,  taking  this  time  into 
consideration,  the  time  for  the  layer  to  reach  concentration 
of  0.2M  is  considered  to  be  53-15  =  3$  seconds.  The  time 
taken  to  allow  2  x  10”  ^  moles  of  1-pentene  to  pass  through 
the  cell  is  only  a  fraction  of  the  time  during  which  the 
layer  has  equilibrium  concentration  of  solute  above  it, 

2/3£  ~  5$;  and  from  this  it  is  concluded  that  the  rate  of 
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depletion  from  the  cell  cavity  is  much  less  than  the  rate 
of  sample  delivery  to  the  cell,  and  that  the  diffusion 
data  obtained  is  valid  for  the  1-pentene — fluorolube 
system.  Because  of  the  similarity  of  vapour  pressures  of 
2-methylbutane  and  diethyl  ether  to  that  of  1-pentene,  these 
conclusions  can  be  extended  to  these  systems.. 

There  was  some  doubt  about  using  acetone  as  a  sample, 
since  it  has  a  low  vapour  pressure  compared  to  the  compounds 
mentioned  above  (100mm  at  7.7°,  400mm  at  39*5°  and  760mm  at 
56.5°)  (4  ).  With  low  vapour  pressures  it  seemed  possible 
that  the  rate  of  delivery  might  be  reduced  and  that  the  rate 
of  depletion  might  become  large  enough  to  affect  the 
diffusion  data.  Calculations  similar  to  those  done  for  the 
l-pentene--f luorolube  system  were  done  for  the  acetone — 
Carbowax  600  system. 

Solutions  of  acetone  in  Carbowax  600  in  known  con¬ 
centrations  were  prepared  (Table  IX) .  The  cell  was 
assembled  with  CC-A  and  the  254  micron  spacer  and  a  layer 
of  the  solution  was  applied.  The  cell  was  inserted  in  the 
spectrophotometer  and  the  wavelength  was  scanned  from  5 
6  microns.  The  peak  height  at  about  5-75  microns  (maximum 
wavelength  for  C=0  stretching  in  ketones)  was  measured  and 
the  resulting  ^A  value  was  calculated  for  each  solution. 

A  Beer’s  law  plot  of  versus  concentration  was  made 
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TaBLE  IX.  ABSORPTION  AT  5. 7 5  MICRONS  OF  SOLUTIONS 
Of  ACETuNc,  IN  CARBOWAX  600  AT  VARIOUS  CONCENTRATIONS 


Solut ion 

Concentrat  Lon 

A  A 

A 

'  0.35  M 

0.065 

B 

0.60  M 

0.12 

C 

0.00  K 

0.15 

Concentration,  M 


Figure  16.  BeerTs  Law  plot,  absorbance  versus  con¬ 
centration  for  solutions  of  acetone  in  carbowax  600 
(Data  frum  Table  IX) 


(fig.  16)  and  the  line  drawn. 


Values  of  vapour  pressures  of  acetone  were  plotted 
at  several  temperatures  (fig.  15)  and  the  value  at  room 
temperature  (  ~  25°c)  was  taken  from  the  curve. 

The  cell  was  assembled  with  CC-A  and  127  micron 
spacer  and  a  layer  of  Carbowax  600  was  aoolied.  With 
acetone  as  sample,  the  water  bath  temperature  at  25°, 
flow  rate  at  0.55  ml/sec,  the  spectrophotometer  wave¬ 
length  selector  locked  at  5-75  microns,  CC-B  was  placed 
on  plate  and  the  procedure  for  measuring  diffusion  co¬ 
efficients  was  followed.  A  point  was  chosen  on  the  re¬ 
sulting  spectrum  at  time  216  seconds  from  zero  time. 

The  -A A  was  calculated,  and  assuming  this  value  repre¬ 
sents  the  average  concentration  of  solute  in  layer,  the 
concentration  was  read  from  the  BeerTs  lav/  plot  (fig.  16). 
From  the  dimensions  of  the  layer  the  amount  of  acetone  in 
the  film  was  calculated  to  be  1.7  x  10“^  moles  (A  -  0.62  , 
c^-3.2  M)  .  From  the  flow  rate  and  vapour  pressure  the 
time  taken  for  1.7  x  10 moles  of  solute  to  pass  through 
the  cell  cavity  was  calculated  to  be  22  seconds.  An 
estimate  of  the  dead-time  was  made  using  a  thin  layer  of 
Carbowax  600  and  found  to  be  about  1$  seconds.  The  time 
then  taken  for  the  layer  to  reach  equilibrium  concentration 
3.21  M  was  approximately  216  -  l£  —  193  seconds:  the  time 
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taken  to  allow  1.7  x  10  ^  moles  of  acetone  to  pass  through 
the  system  was  about. 11$  of  this  time.  From  this  it  was 
concluded  that  for  the  acetone — Carbowax  600  system  a  con¬ 
stant  concentration  of  acetone  is  maintained  above  the 
layer.  This  conclusion  was  assumed  to  cover  the  acetone 
systems  as  well. 

Ill  RESULTS  AND  DISCUSSION 

111  -  1  Experimental  Results 

The  spectra  were  analyzed  by  replotting  the  per  cent 
transmittance  data  as  absorption  data  and  the  best  curves 
vtfere  drawn.  To  each  of  the  curves  tangents  were  drawn  at 
several  points  and  the  slopes  of  these  tangents  were 
measured.  The  logarithms  of  the  slopes  v:ere  taken  and  for 
each  curve  a  plot  of  log  slope  against  time  was  made.  The 
best  straight  line  through  these  points  was  drawn  and  the 
slope  measured.  From  equation  (xv)  it  is  seen  that  the 
slope  of  this  line  is  equal  to  -  'TT'^D/  2.3  h  ,  and  thus 
from  the  slope  the  diffusion  coefficient  can  be  calculated. 

Examples  of  the  original  spectra  (figs.  17,13,  and  19) , 
are  replotted  in  figures  20,  21  and  22  from  transmittance 
values  to  absorbance  values.  On  these  plots  tangents  are 
drawn  and  the  logarithms  of  the  slopes  for  these  examples 
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Figure  17 .  Spectrum  (percent 
transmittance  versus  time)  of 

diethyl  ether  in  fluorolube  at 

(using  CC-B,  254 j\ 

spacer) . 
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Figure  20.  Replot  of  figure  17  absorbance  versus 


time  with  slopes  of  tangents  at  various  times  constructed. 
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Figure  22.  Reolot  of  figure  19  as  absorbance  versus  time  with 
slopes  of  tangents  at  various  times  constructed. 
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TABLE  X.  SLOPES  (dA/dt)  AND  LOGARITHMS  OF  SLOPES  OF  THE 
TANGENTS  TO  THE  ABSORBANCE  VERSUS  TIME  PLOTS* 


Figure 

Time,  sec 

Slope  (dA/dt) 

Log  Slope  (log  dA/dt) k 

20 

35 

0.0066 

-2.18 

40 

0.0040 

-2.40 

45 

0.0032 

-2.50 

50 

0.0013 

-2.74 

55 

0.0010 

-3.00 

60 

0.00066 

-3.13 

21 

150 

0.00063 

-3.20 

130 

O.OOO45 

-3.35 

210 

O.OOO35 

-3.46 

240 

O.OOOI95 

-3.71 

270 

0.00016 

“3.79 

22 

25 

0.0066 

-2.18 

30 

0.0050 

-2.30 

40 

0.0026 

-2.59 

50 

_  O.OOI65 

-2.78 

60 

0.00090 

-3.01 

70 

O.OOO63 

-3.20 

a.  Figures  20,  21,  22. 

b.  Log  (dA/dt)  has  been  shown  previously  to  be  propor¬ 
tional  to  log  (dc/dt)  and  the  two  expressions  are 
used  interchangeably. 
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Figure  23.  Plots  of  logcic/dt  versus  time  for  the 
absorbance  versus  time  plots  shown  in  figures  20,  21, 
22.  The  data  is  taken  from  Table  X. 
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TABLE  XI.  DIFFUSION  COEFFICIENTS  FOR  SYSTEMS  ILLUSTRATED 
BY  SPECTRA3-  AND  CALCULATED  FROM  THE  SLOPES  OF  THE  LOG  dA/dt 

VERSUS  TIME  PLOTSb 


Figure 


-A 


-B 


Corresponding 

System 

ether-f luoro . 
254/a  spacer 

isopentane- 
f luoro lube 

5$4^\spacer 

isopentane- 

fluorolube 

spacer 


—  A  log  dc/dt  / A  t 


.  0.040 


0.005 


0.024 


Diffusion 

Coefficient 

2.4  x  10  ^ 
1.6  x  10  ^ 

3.3  x  10~5 


a.  Figures 

b.  Figure  £3 

•2-1 

c.  Units  for  diffusion  coefficient  are  cm  sec 

d.  Diffusion  coefficients  are  calculated  from  expression 

2  2 

—  ?r  Dt/  2.3  h  which  equals  the  slope  of  plot  log 
(dc/dt)  versus  time. 
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TABLE  Xll.  VALUES  OF  ALog  (dc/dt)/At  AND  DIFFUSION 
COEFFICIENTS  FOR  1-PENTENE  IN  FLUOROLUBE  FOUND 
EXPERIMENTALLY  USING  CC-B  AND  SPACER  THICKNESSES  127  AND 

254  MICRONS 


Diffusion 

Spacer  Thickness  log( dc/dt)  /  A  t  Coefficient  x  10 


0.020 

0.30 

0.009 

0.13 

0.011 

0.17 

0.024 

0.36 

0.027 

0.41 

0.030 

0.45 

0.024 

0.36 

0.030 

0.45 

0.042 

2.5 

0.062 

3.7 

0.045 

2.7 

O.O46 

2.8 

0.020 

1.2 

0.018 

1.1 

0.022 

1.3 

0.030 

1.3 

0.044 

2.7 

0.043 

1.7 

a 


Units  cm  /  sec 
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TABLE  Xlll.  VALUES  OF  ALog(dc/dt)/  At  AND  DIFFUSION 


COEFFICIENTS  FOR  1-PENTENE  IN  FLUOROLUBE  FOUND 
EXPERIMENTALLY  USING  CC-B  AND  SPACER  THICKNESSES  38I 

584  MICRONS 

AND 

Spacer  Thickness 

"  Alog(dc/dt)/  At 

Diffusion 

Coefficient  x 

105 

381 AA 

0.024 

3.3 

0.012 

1 . 6 

0.027 

3.6 

0.027 

3.6 

0.024 

3.3 

0.008 

1.1 

0.016 

2.2 

0.033 

4.3 

564  yx 

0.005 

1.6 

0.011 

3.5 

0.009 

2.9 

0.012 

3. a 

0.012  • 

3. a 

0.014 

4.5 

0.012 

3. a 

0.007 

2.2 

• 

0.005 

1.6 

0.012 

3. a 

2 

a.  Units  cm  /sec 


. 


-66- 


TABLE  XIV.  VALUES  OF  JS  Log(  dc/dt)  /  A  t  AND  DIFFUSION 
COEFFICIENTS  FOR  2  METHYLBUTANE  IN  FLUOROLUBE  FOUND 
EXPERIMENTALLY  USING  CC-B  AND  SPACER  THICKNESSES  127  AND 

254  MICRONS 


Spacer  Thickness 


254  JA 


Diffusion 

~~  A  log(  dc/dt ) /  At  Coefficient  x  10^ 


0.047 

0.54 

0.029 

0 . 44 

0.022 

0.34 

0.054 

0.31 

0.032 

0.43 

0.023 

0.35 

0.030 

0.45 

0.014 

0.22 

0.020 

0.30 

0.034 

0.51 

0.042 

2.5 

0.023 

1.7 

0.030 

1.3 

0.034 

2.1 

0.023 

1.7 

0.040 

2.4 

0.023 

1.7 

0.033 

2.3 

0.027 

1.6 

a 


Units  cm/ sec 
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TABLE  XV.  VALUES  OF  &Log(dc/dt)/  At  AND  DIFFUSION 
COEFFICIENTS  FOR  2 -ME TH Y LBU T AN E  IN  FLUOROLUBE  FOUND 


EXPERIMENTALLY  USING 


CC-B  AND  SPACER  THICKNESSES  3^1 
534  MICRONS 


AND 


Diffusion 


Spacer  Thickness  — ~  A  Iog( dc/dt ) /  At  Coefficient  x  10 • 


0.024 

3.3 

0.014 

1.9 

0.024 

3.3 

0.026 

3.5 

0.023 

3.1 

0.020 

2.7 

0.019 

2.6 

0.016 

2.2 

0.006 

1.9 

0.007 

2.2 

0.004 

1.3 

0.005 

1.6 

0.004 

1.3 

0.003 

1.0 

a 


2 

Units  cm  /sec 
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TABLE  XVI.  VALUES  OF  A  Log ( dc/dt ) / A t  AND  DIFFUSION 
COEFFICIENTS  FOR  DIETHYL  ETHER  IN  FLUOROLUBE  FOUND 
EXPERIMENTALLY  USING  CC-B  AND  SPACER  THICKNESSES  127  AND 

254  MICRONS 


Spacer  Thickness 


~~ A  log ( dc/dt)  /  A  t 

0.036 
0.022 
0.026 
0 .02$ 

0.034 

0.024 

0.034 

0.017 

0.017 

0.030- 

0.022 

0.039 

0.032 

0.040 

0.034 

0.027 

0.036 

0.040 

0.032 


Diffusion 
Coefficient  x  10^ 

0.54 

0.33 

0.39 

0.42 

0.51 

O.36 

0.51 

0.26 

0.26 

0.45 


1.3 

2.4 
1.9 

2.4 
2.1 
1.6 
2.2 

2.4 
1.9 


a 


Units  cm  7  sec 
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TABLE  XVII.  VALUES  OF  A  Log  (dc/dt)/A  t  AND  DIFFUSION 
COEFFICIENTS  FOR  DIETHYL  ETHER  IN  FLUOROLUBE  FOUND 
EXPERIMENTALLY  USING  CC-B  AND  SPACER  THICKNESSES  331  AND 

564  MICRONS 


Diffusion  a 

5 

Spacer  Thickness  —A log  (dc/ dt)/ A  t  Coefficient  x  10 


381 

0.013 

1.8 

• 

0.017 

2.3 

0.019 

2.6 

0.023 

3.1 

0.014 

1.9 

0.013 

2.4 

0.017 

2.3 

584 

0.011 

3.5 

0.011 

3.5 

- - 

0.006 

1.9 

-  — 

0.007 

2.2 

0.009 

2.9 

0.010 

3.2 

a .  Units  cm  / sec 
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TABLE  XV111.  SUMMARY  OF  VALUES  OF  DIFFUSION  COEFFICIENTS 
FOR  1- PEN TEN E ,  2 -ME TH YLBU T AN E  AND  DIETHYL  ETHER  IN  FLUOROLUBE 


Solute 

Average  Diffusion 

Coefficient  x  10^a 

Average 

Deviation  x  10 ^ 

Standard 

5 

Deviation  x  10 

1-pentene 

2.7  - 

0.9 

1.0 

2-Me-butane 

2.2 

0.6 

0.7 

di-Et-ether 

2.4 

0.4 

0.6 

Units  cm  y sec 
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TABLE  XIX.  VALUES  OF  DIFFUSION  COEFFICIENTS  FOR  ACETONE 
CARBON  AX  600,  CARBON' AX  1500,  SILICONE  GREASE,  DINONYL- 


PHTHALATE 

AND  UCON 

Spacer 

Diffusion 

Solvent 

Thickness 

5  2 

Coefficient  x  10  cm  /sec 

Carbowax 

600 

127yu 

0.12 

0.060 

0.070 

381 r 

1.5 

2.1 

1.2 

Carbowax 

1500 

127ia 

0.030 

0.040 

Silicone 

Grease 

127  yx 

0.12 

0.060 

% 

0.98 

Dinonyl- 

127yx 

0.068 

phthalate 

0.15 

0.18 

r-f" 

to 

1.7 

1.7 

Ucon  (-50-HB 
280X) 

127 

0.065 

0.060 

0.15 

38lyU. 

1.9 

1.7 

•  * 
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are  presented  in  lable  X.  The  plots  of  log  slope  versus 
time  are  shown  in  figure  23  and  the  results,  i.e.,  slopes 
of  lines  and  calculated  diffusion  coefficients,  are  shown 
in  Table  XI. 

The  slopes  of  the  log  dc/dt  versus  t  plots  and  the 
calculated  diffusion  coefficients  for  the  1-pentene-- 
fluorolube,  2-methylbutane — fluorolube  and  diethyl  ether — 
fluorolube  systems  are  presented  in  Tables  Xll  -  XV111 
inclusive.  Final  tabulation  of  results  as  averages, 
average  deviations  and  standard  deviations  for  the  fluoro¬ 
lube  systems  is  given  in  Table  XV111. 

From  the  spectra  it  was  seen  that  as  the  spacer 
thickness  increased,  the  change  in  absorption  or  the 
equilibrium  concentration  decreased. 

'  Similar  data  and  results  for  the  acetone — Carbowax 
600,  Carbowax  1500, _ silicone  grease,  DNP  and  Ucon  systems 
are  presented  in  Table  XIX. 

111-2  Discussion 

The  theoretical  model  of  diffusion  into  a  slab  coupled 
with  Fickrs  second  lav/  of  diffusion  lead  to  a  simple 
expression  for  the  diffusion  coefficient  as  a  function  of 
the  change  of  concentration  with  time  at  a  given  distance 
in  the  slab.  This  change  of  concentration  is  easily 
measured  by  internal  reflectance  spectroscopy.  So  it  appears 
that  the  method  outlined  in  this  work  should  yield  fairly 
accurate  values  for  the  diffusion  coefficients.  However, 


'll 


the  absence  of  diffusion  data  in  the  literature,  especially 
for  the  systems  such  as  those  present  in  gas  chromatography 
makes  it  difficult  to  affirm  the  values  of  D  derived  experi 
mentally . 

The  data  which  appears  in  Tables  X11-XV11  inclusive 
shows  the  discrepancy  in  values  measured  for  all  samples 
using  the  127  micron  spacer  with  respect  to  those  values 
measured  at  higher  spacer  v/idths .  The  data,  which  showed 
lower  values  for  the  diffusion  coefficients,  was  deleted 
for  the  final  analysis  and  was  not  entered  into  Table 
XV111.  There  are  two  possible  sources  of  error  for  thin 
layers,  i.e.,  non-uniformity  of  the  film  and  surface  effect 
The  uniformity  of  the  film  cannot  be  guaranteed,  par¬ 
ticularly  not  for  very  thin  layers  of  liquids.  Liquids 
such  as  fluorolube,  DNP  and  Ucon  can  flow,  and  whether 
this  flow  is  caused  by  movement  of  the  cell,  or  by  the 
flow  of  gas  in  the  cavity,  this  flow  of  liquid  is  more 
important  for  thin  layers.  If  the  solute  lowers  the 
surface  tension  of  the  solvent,  the  solute  will  tend  to 
be  more  concentrated  at  the  surface  of  the  layer.  This 
effect  would  give  values  of  the  diffusion  coefficient 
lower  than,  expected.  Since  thicker  layers  appeared  to 
be  more  stable  in  preparation  and  since  the  surface 
effect  was  expected  to  have  smaller  influence  on  the 
results  for  thicker  layers,  only  data  obtained  from 
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measurements  with  254^,  3#1 and  5$4yufrom  Tables  Xll- 
XVII  were  considered  in  the  calculation,  of  diffusion  co¬ 
efficients  . 

The  author  could  not  explain  the  aprearance  of  the 
spectra  with  regard  to  decreasing  equilibrium  concentra¬ 
tion  with  increasing  spacer  thickness  as  mentioned  in 
section  111-1#  A  study  was  made  using  a  solution  of  known 
concentration  to  see  if  spacer  thickness  had  any  effect 
on  the  absorption.  With  all  spacers  tried  the  absorption 
was  the  same.  This  is  also  expected  in  the  measurement 
of  diffusion  coefficients,  i.e.,  the  equilibrium  concen¬ 
tration  should  be  independent  of  spacer  thickness.  Thus 
no  explanation  f.or  this  anomalous  behaviour  is  offered. 

The  results  for  each  solute-f luorolube  system  v/ere 
averaged  and  the  average  deviation  and  standard  deviation 
calculated.  These  results  were  recorded  in  Table  XV111. 
Because  the  standard  deviation  is  large,  the  diffusion 
coefficients  measured  cannot  be  considered  quantitative. 
However  some  qualitative  observations  can  be  made. 

The  Wilke-Chang  equation,  considered  by  Giddlngs  (10) 
to  be  one  of  the  best  empirical  expressions  for  diffusion 
coefficients,  D=  7.4  x  10  ^  (xM)2T/yv^*^  was  used  to 
try  to  estimate  the  predicted  values  for  the  1-pentene — 
fluorolube  system.  The  viscosity  of  fluorolube  at  100°F 
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[  3$  C)  (27)  was  100  centistokes  -  as  centistokes  times 

density  equals  centipoise,  and  the  density  is  given  for 
100°F  as  1.92  gm/ml,  the  viscosity  in  centipoise  is  found 
to  be  192.  The  average  molecular  weight  of  fluorolube  was 
taken  to  be  775*  Taking  the  molal  volume  of  1-pentene  as 
111  (from  values  in  Table  IV) ,  the  Wilke- Chang  equation 
gives  the  diffusion  coefficient  at . 100°F  as  about  2  x  10 ”7 

cm2/Sec.  The  value  at  25°C  will  be  less  than  2  x  10 
because  of  the  increase  in  viscosity  with  decrease  in 
temperature.  The  Wilke-Chang  estimate  is  lower  than  the 
experimental  values  found.  A  similar  case  was  reported 
in  the  literature  when  the  values  for  the  diffusion  co¬ 
efficient  for  the  isobutylene--DNP  system  investigated 
by  Houghton  et.  al.  (15)  ,  D  for  temperature  range 

0  -  75°C  were  0.026  -  0.61  x  KPcm^/sec,  we re  higher  than 
the  Wilke-Chang  estimates  by  factors  of  11  to  44*  Giddings 
writes  (12)  "Thus  until  further  evidence  is  available, 
calculations  regarding  small  solute  molecules  in  large 
molecules  of  gas  liquid  chromatography  are  suspect." 

Although  the  Wilke-Chang  equation  does  not  give  a 
good  estimate  of  the  experimentally  measured  diffusion  co- 

r 

efficients,  qualitative  relationships  can  be  made  using 
the  Wilke-Chang  proportionality  D  oc  l/V0-6.  The  molal 
volumes  for  2-methylbutane  and  diethyl  ether  were  calcu¬ 
lated  from  Table  IV  (113  and  107  respectively),  and  the 
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values  of  experimental  diffusion  coefficients  and  range 
defined  by  the  average  deviations  v/ere  plotted  against  the 
values  of  l/V  (fig.  24) .  Although  the  three  points  do 
not  form  a  straight  line,  considering  the  range  of  values 
caused  by  the  deviation  there  appears  to  be  a  trend  to 
show  that  the  diffusion  coefficient  is  a  linear  function 
1/V^  *  ^  and  that  if  a  large  number  of  solutes  of  varying 
molal  volume  were  run  in  fluorolube,  a  variation  of  the 
Wilke-Chang  equation  might  be  obtained. 

Of  the  acetone-solvent  series  (solvents  -  Carbowax 
600,  Carbowax  1500,  DNP,  etc.)  the  values  of  diffusion 
coefficients  were  made  from  measurements  with  both  the 
127  and  3#1  micron  spacers  for  the  solvents  of  lower 
viscosity  (Carbowax  600,  DNP,  Ucon) ;  and  for  the  solvents 
of  higher  viscosity  (Carbowax  15.00,  Silicone  grease)  runs 
were  only  made  with  the  127  micron  spacer.  Only  a  few 
runs  were  made  and  thus  no  averages  were  taken. 

Since  viscosity  and  molecular  weight  data  were  not 
available  for  all  solvents  -  only  available  for  the 
acetone — DNP  system  -  and  since  the  Wilke-Chang  equation 
was  shown  to  be  unreliable  from  a  quantitative  consider¬ 
ation,  only  qualitative  comparisons  could  be  maae.  Using 
the  data  for  acetone — DNP,  acetone  Carbowax  600,  and 
acetone — Carbowax  1500,  comparison  was  mode  on  the  basis 
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of  the  Wilke-Chang  dependence  of  D  on  Ms/h  and  estimates 
of  the  viscosities  of  the  carbowaxes  were  made.  The 
diffusion  coefficient  for  the  acetone — DNP  system  was 
found  to  be  1.7  x  10  5  cm^/sec  and  was  about  the  same  as 
the  value  for  diffusion  of  acetone  in  Carbowax  600 
(average  1.6  x  10  5  cm sec) .  The  value  for  this 
latter  system  is  about  fifty  times  as  large  as  the  diffu¬ 
sion  coefficient  for  the  Carbowax  1500  system.  From  the 
viscosity  of  dinonylphthalate ,  Si  centipoise  at  25°C  (15), 
and  the  molecular  weight  the  ratio  M2/0  was  found  to  be 
about  0.25.  Thus  this  ratio  for  the  Carbowax  600  system 
should  be  the  same  and  for  the  Carbowax  1500  system  it 
should  be  0.005  -  if  the  Wilke-Chang  dependancy  is  appli¬ 
cable.  From  this  ratio  the  viscosities  of  the  carbowaxes 
were  estimated:  for  Carbowax  600  viscosity  ^100,  for 
Carbowax  1500  viscosity  4700  (values  in  centipoise) . 
Although  the  viscosities  for  the  carbowaxes  are  not  known 
at  25°C,  the  values  found  by  using  the  Wilke-Chang  de¬ 
pendence  of  diffusion  coefficient  on  molecular  weight  and 
viscosity  give  values  which  appear  to  be  in  agreement  with 
the  properties  observed  while  handling  these  reagents. 
Carbowax  600  is  liquid  with  pouring  qualities  much  like 
those  of  DNP  as  is  expected  with  viscosity  of  100  centi¬ 
poise  ;  Carbowax  1500  is  a  waxy  solid  which  is  also  expected 
from  its  high  value  of  viscosity. 
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In  the  measurement  of  diffusion  coefficients  and 
qualitative  comparisons  with  the  YVilke-Chang  equation  no 
mention  was  made  of  temperature  effects.  The  dependence 
of  diffusion  coefficients  on  temperature  in  the  YJilke- 
Chang  relation  is  shown  to  be  two  fold,  i.e.,  D  is 
directly  proportional  to  T,  and  the  D  is  inversely  pro¬ 
portional  to  viscosity  which  is  in  turn  inversely  pro¬ 
portional  to  the  temperature.  The  change  of  viscosity 
with  temperature  is  great  especially  for  high  molecular 
weight  solvents.  Thus  any  change  in  temperature  during 
the  course  of  the  experiment  would  cause  a  change  in  the 
value  of  D  calculated.  Since  with  infrared  spectroscopy 
the  sample  heats  up  during  the  scan  (more  for  thick 
layers  than  for  thin) ,  the  values  of  the  diffusion  co¬ 
efficients  would  be  higher  than  those  expected  from  the 
Wilke-Chang  equation. 

These  comparisons  offer  a  qualitative  confirmation 
of  the  values  of  diffusion  coefficients  found  in  this  work. 

111-3  Errors 

There  are  sources  of  error  which  are  inherent  in  the 
method ; ✓ the  replotting  of  data  of  percent  transmittance  as 
absorption  data,  the  construction  of  curves  and  tangents, 
and  the  measurements  of  slopes.  The  error  in  replotting 
data  is  dependent  on  the  value  of  T  being  converted. 
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The  spectra  could  be  read  to  the  nearest  0 . 5^T  and  in  Table 
XX  percent  error  in  absorbance  is  shown  at  various  values  of 
$T .  From  this  table  it  can  been  seen  that  at  high  and  low 
values  of  yoT  the  percent  error  in  absorbance  rises.  The 
error  in  the  construction  of  curves  and  tangents  and  in 
measuring  slopes  was  estimated  by  having  the  data  re¬ 
plotted  independently  and  comparing  the  slopes  of  the  log 
(dc/dt)  versus  time  plots  with  those  measured  by  the  author. 
The  results  of  this  replotting  and  comparison  with  the 
authorTs  results  are  shown  in  Table  XXI  along  with  an  estimate 
of  the  percent  error  between  the  two  sets  of  results. 

The  main  source  of  error  in  this  procedure  arises 
from  the  preparation  of  the  solvent  layer.  There  can  be 
no  guarantee  as  to  the  uniformity  and  actual  thickness  of 
the  layers,  especially  when  using  the  thinner  spacers. 
Overflow,  and  leakage  under  the  spacer  due  to  non-uniformity 
of  the  grease  seal  can  be  controlled  to  some  extent,  but 
these  can  occur  causing  variation  in  the  layer  thickness 
and  consequently  variation  in  the  value  of  the  diffusion 
coefficient.  N 

Another  source  of  error  which  was  mentioned  previously 
is  temperature  and  the  heating  of  the  solvent  layer  during 
the  spectroscopic  measurement  of  diffusion  coefficients. 

One  way  in  which  this  could  be  controlled  would  be  the  use 
of  a  filter  which  would  cut  out  the  radiation  less  than 


TABLE  XX.  ERRORS  IN  ABSORBANCE  DUE  TO  TRANSFORMATION 
OF  DATA  FROM  PERCENT  TRANSMITTANCE  TO  ABSORBANCE  UNITS® 


ins  % T 

Absorbance 

Absolute  error 

Percent  error 

5 

1.301 

0.040 

3 . 1  io 

10 

1.000 

0.023 

■  2.4/ 

20 

0.699 

0.009 

1  •  3/ 

30 

0.522  • 

0.006 

1.2/a 

50 

0.301 

0.004 

1.4/ 

70 

0.155 

0.003 

2.0 / 

90 

0.045 

0.004 

9.7/ 

a  A  -  I03  t/T 


two  microns  (the  bulk  of  the  energy  emitted  by  the  black 
body  source. 

Despite  these  sources  of  error,  this  method  does  give 
reasonable  values  of  diffusion  coefficients  ana  appears 
to  be  a  good  method  for  determining  them.  The  accuracy  of 
these  numbers  could  be  improved  by  changing  aspects  of  the 
procedure  and  cell  design.  The  use  of  thicker  spacers  - 
of  the  order  of  one  millimeter  -  should  reduce  many 
problems  in  non-uniformity  of  layer  and  possible  reduction 
of  surface  effects.  These  thicker  spacers  should  also 
make  it  easier  to  apply  a  more  uniform  and  reproducible 
layer  of  solvent.  The  use  of  teflon  as  spacer  material 
instead  of  metal  would  improve  the  contact  between  the 
prism-plate  and  the  spacer  and  would  eliminate  the  need 
for  a  grease  seal.  If  the  spacer  were  glued  onto  the 
plate,  the  need  for  a  cell  cover  such  as  CC-A  would  be 
eliminated  and  the  length  of  time  required  for  an  indi¬ 
vidual  run  would  be  decreased. 
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TABLE  XXI. 


ERRORS  FOUND  IN  REPLOTTING  OF  DATA 


Slope  of  log  (dc/dt) 

versus  t  plot  Percent  Error 


Measurement 

1 

11 

A 

0.007 

0.008 

14 

B 

0.034 

0.028 

21 

C 

0.011 

0.0085 

34 

D 

0.028 

0.024 

17 

E 

0.014 

0.017 

21 

F 

0.014 

0.017 

21 

G 

0.026 

0.023 

13 

H 

0.027 

0.023 

17 

I 

0.028 

0.032 

14 
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IV  CONCLUSION 

Using  theoretical  background  of  diffusion  as  given 
by  Fick!s  law  and  dost  (16),  and  the  technique  of  FMIR, 
a  method  for  the  determination  of  liquid  diffusion  co¬ 
efficients  was  worked  out.  Two  cell  covers  for  the  FMIR 
solid  sampling  cell,  CC-A  and  CC-B,  were  designed  to  pro¬ 
vide  uniform  liquid  films  of  known  thicknesses;  and  the 
sampling  system  was  designed  to  provide  constant  flow  of 
sample  in  the  gas  phase  to  the  layer.  Using  a  Perkin- 
Elmer  IR-137  on  its  side  to  allow  horizontal  placement 
of  the  prism-plates,  diffusion  coefficients  for  systems 
of  different  solutes  in  single  solvent  and  one  solute  in 
different  solvents  were  measured  by  following  the  change 
of  concentration  with  time  as  indicated  by  the  spectra. 
Although  the  data  cannot  be  compared  to  literature  values, 
they  are  consistent  with  approximations  given  by  authors 
such  as  Perrett  (19)  and  Bohemen  and  Purnell  (3)*  The 
values  obtained  and  that  could  be  obtained  by  this  method 
will  help  elucidate  liquid  phase  mass  transfer  terms  in 
gas  chromatography. 
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